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Abstract. We report the statistical properties of classical particles 
in (2-1-1) gravity as resulting from numerical simulations. Only particle 
momenta have been taken into account. In the range of total momentum 
where thermal equilibrium is reached, the distribution function and the 
corresponding Boltzmann entropy are computed. In the presence of large 
gravity effects, different extensions of the temperature turn out to be in- 
equivalent, the distribution function has a power law high-energy tail and 
the entropy as a function of the internal energy presents a flex. When the 
energy approaches the open universe limit, the entropy and the mean value 
of the particle kinetic energy seem to diverge. 



1. Introduction. The inclusion of gravitational effects in thermodynamics has 
important consequences at the cosmological level. There is a particular interest in this 
subject also because gravity deeply modifies certain basic features of non-gravitating 
systems; for instance, in the presence of gravitational interactions one cannot define 
the thermodynamic limit in the standard way, one may find bounds for the value of the 
entropy, the additivity of the thermodynamic potentials is no longer valid in general, 
etc. In this article we explore the statistical properties of a classical gravitating gas in 
(2+1) dimensions with vanishing cosmological constant. The gravitational dynamics 
in (2+1) dimensions is rather simple but it is not trivial, so this is one of the simplest 
models in which one can look for new phenomena. 

Our results are based on numerical simulations. We have considered a set of clas- 
sical pointlike particles, with random initial momenta, moving on a spatial surface of 
topological type either or the sphere for a finite time interval. We followed only 
the evolution of the particle momenta. Let G be the gravitational constant and U the 
total energy of the system. When AGU < 1 , the system reaches thermal equihbrium. 
Whereas, for space-hke total momentum or when AGU > 1 , the numerical analysis 
shows the presence of dynamical instabilities; jets of particles appear with no-limits 
increasing energy. 

In the case AGU < 1 , the resulting distribution function in momentum space / 
is numerically produced and the corresponding Boltzmann entropy, yS" oc J f log / , 
is computed as a function of the total energy U of the system. If AGU <^ 1 the 
gravitational effects can be neglected and the entropy S shows the classical logU 
behaviour. On the other hand, when the value of the energy approaches the limit 
AGU = 1 , gravity has a strong influence on the dynamics of the particles and the value 
of S seems to diverge. The implications of this behaviour on the heat capacity of the 
gas are considered. Finally, we compare various definitions of temperature which turn 
out to be inequivalent in the energy range of large gravity effects. 

2. Classical gravitating particles. The peculiar properties of (2+1) gravity have 
been described in Ref. [l,-,4]. In any (2+1) spacetime M containing classical pointlike 
particles, all the curvature is located on the worldlines of the particles. Each parti- 
cle moves in a locally fiat surrounding spacetime. Since the isometry group of flat 
spacetime is isomorphic with the Poincare group IS0{2, 1), both the kinematics and 
the laws of the dynamics are covariant with respect to transformations of the Lorentz 
group -50(2,1). 

A rather simple description of the spacetimes containing N pointlike gravitating 
particles has been produced by 't Hooft [5,6]. This description is given by the lin- 
ear evolution of Cauchy surfaces which are tiled by spatial planar polygons; the time 
evolution also includes modifications of the tiling combinatorics according to codified 
transition rules [6]. A Lorentz frame is associated to each polygon with local coor- 
dinates {t, X, y) and metric ds^ — dt^ — dx^ — dy^ ; with respect to this frame, the 
interior of each polygon belongs to the spatial t — constant surface. The extrinsic 
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curvature is vanishing in the interior of each tile and is singular on the edges. Each 
particle is placed in one corner of a polygon; this corner is the intersection point of two 
consecutive edges which must be identified. 

At a given initial time t = to , one can always introduce a tessellation in which 
all the particles are placed at the corners of one polygon V{to) ; let us denote by LS 
the Lorentz system which is associated with V{to) . The initial values {p^-^ } (with 
i = 1,2, ■■■,N) of the particle momenta are defined with respect to the frame LS . In 
addition to V{to) , the (t = to) tessellation may possibly contain several polygons; the 
complete description of the initial configuration includes all the lengths and velocities 
of the polygon edges. 

During a short time interval after t — to , one has a "linear" evolution. Because 
of the constant velocities of the edges, the length of each edge is a linear function of 
time, the angles between edges and the particle momenta are constant in time. When 
the length of one edge vanishes or when one of the corners hits an edge, one has a 
transition in which the structure of the tessellation undergoes a local change. So, the 
nontrivial part of the evolution consists of a sequence of transitions; in each transition 
the shape and the number of polygons, together with the velocities of the edges, are in 
general modified according to a given set of rules [6]. 

After a finite time interval At = ti — to in which no big-crunch occurs, one finds 
a tessellation in which the particles are in general placed at the corners of different 
polygons. In order to compare the initial {t — to) and the final {t — ti) values of 
the particle momenta, we need to refer these values to a unique coordinate system; we 
shall choose the Lorentz frame LS as reference system. 

Let Qo be a point in the interior of V{to) ; we shall denote by Qi the point in M 
which has the same spatial coordinates as Qq and has time coordinate ti with respect 
to the frame LS . We assume that Qi does not belong to the worldline of a particle; 
for, if this is the case, we can simply modify the spatial position of Qq so that the 
corresponding Qi does not belong to a worldline. In a neighbourhood Q G M of Qi 
the spacctimc is flat and one can introduce local coordinates which arc adapted to the 
Lorentz system LS . The piece in fl of the spatial surface t = ti can be extended 
in M and one will eventually reach the worldlines of the particles. This extension is 
not unique, in general. In order to remove all ambiguities, we shall require that the 
extended surface t = ti is star-shaped with respect to the geodesic lines of minimal 
(spatial) distance emerging from Qi ; this simply means that if a point Q belongs to 
the surface, then all the points of the (minimal-length) geodesic connecting Q with 
Qi also belong to the surface. This extension of the spatial surface t — ti determines 
a polygon V{t-i) with the property [6,7] that each particle is placed at a corner of 
V{ti) . Moreover, one can introduce a new (t = ti) tessellation of space containing the 
tile V{ti) . With respect to the Lorentz frame LS of V{ti) , which is the same frame 
of V{to) , the final {t = ti) values of the particle momenta will be denoted by { k'^-^ } . 

To sum up, with respect to a fixed Lorentz system, the time evolution results in a 
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modification of the particle momenta 

U{h,to) : {pl)} {k^)} , 

where the map f/(ti,to) depends on the initial {t = to) data. We shall conclude this 
section by recalling the general structure of U{ti,to) ■ 

Consider a set of oriented closed paths { 7i } on V{to) which are based on the 
point Qq and have no intersection with the particle worldlines. Let the loop 7^ 
be simply linked with the worldhne of the i-th particle; the paths { 7? } are closed 
because some of the edges of 'P(to) must be identified in the tessellation. The effects 
of a parallel transport of a vector along the loop 7,4 can be described by means of a 
Lorentz transformation i?(7i) acting on the tangent space in Qq . Since M is fiat in 
a neighbourhood of Qo , the Lorentz frame LS canonically defines a reference system 
in the tangent space at Qq . The spacetime curvature is concentrated on the particle 
worldlines, so H{^i) is invariant under smooth deformations of 7j . Moreover, one has 
-^(7i~^) — H'^i^i) where 7"^ represents 7j with reversed orientation. Now, with 
a suitable choice of the loops { 7i } and a given (clockwise) prescription for the loop 
orientations, one finds [6] 

where { Ja } ^-re the generators of the Lorentz group and { p"^^ } are the components of 
the momentum of the i-th particle with respect to the reference system LS . Since all 
the particles are placed at the corners of the polygon P(to) whose boundary -without 
edges identifications- has the topology of , an ordering of the particles turns out to 
be fixed modulo cyclic permutations; let us assume that the natural ordering of the 
index i = 1,2, ...,N corresponds to consecutive particles on the boundary of V{to) ■ 
The initial (t = to) particle data can conveniently be expressed by means of an ordered 
set / of SO {2, 1) group elements 

I = {//(7i),i/(72),---,/^(7iv)} = {e^w,e?(^),...,e^(-)} . 

We now need to introduce the Lancaster-Sasakura representation [8] of the braid group 
Bn acting on I ; the generator aj oi (with j — 1,2, N — 1) is represented by 

aj : { e^w , . . . , e^(^) } ^ { e^''^ , • • • , e^^^"'' , e^^^+'^ , e-^o+^e^^t^'e^o+i) , • • • , e^(^) } . 

Let us denote by F the ordered set of SO (2,1) group elements 

F = { (h) , e^(2) ^ . . . ^ efc(iv) j 

where { /c^^^ } are the final {t — ti) particle momenta. One can show [6,-, 9] that the 
time-evolution map is given by an element g{ti,to) e in the Lancaster-Sasakura 
representation 

U{ti,to) : I ^ F ^ 9{h,to)I . (2.1) 
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The particular element g{ti,to) e Bjsr which enters equation (2.1) depends on the 
initial positions and velocities of the particles. Any element of the braid group can be 
written as an ordered product of generators. So the time evolution can be interpreted 
as the result of a certain number of consecutive elementary processes; each elementary 
process is associated with a generator of the braid group. 

3. Computation rules. For any specific choice of the initial lengths and edge ve- 
locities of the starting tessellation, 't Hooft rules determine a unique time evolution. 
With a large number of particles, the complete evaluation of the exact microscopic 
state of the system is rather laborious [6]. On the other hand, we are interested in 
the thermodynamic properties of a large number N of gravitating particles. There- 
fore we shall introduce certain approximations which simplify the computation of the 
macroscopic variables. The main idea consists of adopting a iterative method -based on 
some random process- to determine the Braid group element g{ti,to) which specifies 
the evolution (2.1) of the particle momenta. 

If one has a 3-dimensional classical gas in standard conditions of temperature and 
pressure, one can approximate the molecular motion by a sequence of diffusion processes 
among randomly chosen pairs of molecules. This is not an accurate approximation at 
the microscopic level but turns out to be a rather good method to evaluate macroscopic 
variables. Similarly, for the 2-dimensional gravitational gas we will approximate the 
time evolution by a sequence of "scattering processes" involving randomly chosen pairs 
of particles. 

Suppose that a and b are two integers which have been randomly chosen among 
the set {1,2, ...,iV} with the condition a ^ b. Suppose that, for instance, a < b. 
One can then define a particular element h{a,b) of the braid group which is the 
product of {b — a) generators 

h{a, b) = ab-i ■ ■ ■ aa+i ■ (Ja 

The elementary process associated with the integers a and b can be interpreted as 
a scattering process involving the a-th and 6-th particles. After this scattering, the 
particle momenta are described by the new ordered set of 5*0(2, 1) elements 

I' = h{a,b)I , (3.1) 

or 

I' = h~\a,b)I . (3.2) 

Possibilities (3.1) and (3.2) are assumed to have equal probabilities to occur. After this 
first step, one has simply to repeat the procedure: one determines two new random 
numbers a' and b' and applies h{a',b') (or h~^{a',b'), with probability 1/2) to the 
ordered set /' in order to find /" , and so on. 

One can use a simpler recipe to approximate the time evolution. At each step of 
this new recipe, one randomly determines an integer a e {1,2,...,A^ — 1} and the 
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modification of the particle momenta is described, with equal probabilities, by the 
generator aa or its inverse in the Lancaster-Sasakura representation. We have 
verified that, in all the considered examples, both methods lead to the same conclusions. 

Our computation rules only concern the values of the particle momenta; the spa- 
tial positions of the particles and their mutual distances have been ignored. In the 
ordinary 3-dimensional case of a gas inside a finite volume, one can safely ignore the 
spatial positions of the molecules because the distribution of equilibrium is in fact ho- 
mogeneous. So, as long as the gas is approaching the thermal equilibrium, the actual 
position of each molecule plays no role in the evaluation of the macroscopic variables. 
For this reason, the assumption of a homogeneous equilibrium distribution could be 
introduced to justify our computation rules. Actually, the 2-dimensional gravitating 
gas that we consider is not evolving inside a container with fixed walls; consequently, 
the equilibrium particle distribution -when it exists- is not necessarily homogeneous 
in space. Yet, the computation rules that we use to approximate the evolution of the 
particle momenta are expected to be valid anyway because each elementary process in- 
volving two gravitating particles in (2-1-1) dimensions does not depend on the distance 
between these particles. 

4. Conservation law. Since we ignore the spatial positions of the particles as well 
as the whole "extension" of the system, we can only compute how the thermodynamic 
potentials depend, for instance, on the internal energy U of the system. So we need 
to give a definition of U and discuss the associated conservation law. 

The components = {p^,p) of the momentum of a classical pointlike particle of 
mass m can be determined by measuring the velocity of the particle with respect to 
a local Lorentz frame. Because of Lorentz covariance, the energy p° and the spatial 
components p of the momentum have the usual dependence on the mass and on the 
velocity of the particle. One has (p)^ = p"Pa = m'^ ■ For a system of A^ > 1 gravitating 
particles, the sum of the energies of the particles is not a conserved quantity. In 
the presence of gravitational interactions, the ordinary conservation law for the 3- 
momentum gets modified; the new conservation law [7] for the 3-momentum is Lorentz 
covariant and can be expressed in terms of 5*0(2, 1) group elements. Since the time 
evolution of the particle momenta is given in equation (2.1), it is easy to verify that, 
whatever g{ti,to) € Bn might be, the composite element 

e?(i) . • . • e^w = e^^*"*) = exp ( SttG p^^^^) J„ ) (4.1) 

is conserved; { p^^^j) } are the components of the total 3-momentum of the A^ particles. 
In (2-1-1) dimensions the gravitational field does not describe propagating degrees of 
freedom; indeed, the total 3-momentum only depends on the variables associated with 
the particles. No Newtonian or potential energy is present in (2+1) dimensions, and 
in fact { p^jo^) } do not depend on the positions of the particles. As expected, in 
the G — > limit, expression (4.1) reproduces the usual additive conservation law of 
3-momentum which is valid in ordinary Minkowski space. 
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The total 3-momentum admits an intrinsic definition. Consider an oriented loop 7 
in M based on the fixed point Qq & M and let E be an oriented disc in M whose 
boundary coincides with 7 . Each particle worldline has a natural orientation and its 
intersection points with S have definite signatures. We will say that a given worldline 
is enclosed by 7 when the algebraic sum of its intersections with E is equal to unity. 
Suppose now that all the worldlines of the particles are enclosed by 7 ; the path-ordered 
integral of the spin connection (which must be defined in a smooth coordinate system) 
along 7 determines the total 3-momcntum. In facts, the 5*0(2, 1) element which 
describes the effects of a parallel transport of vectors along 7 is given by 

H{^) = exp(87rGpf,„,)J„) . (4.2) 

When the total 3-momentum is timelike, one can find a Lorentz transformation A such 
that 

A • exp ( SnG p^^^j) J„ ) • A"^ = exp ( SnG UJq ) , (4.3) 

where U is called the internal energy of the particles system. The generator of rotations 
Jo has integer eigenvalues and then 8nGU is apparently defined only modulo 27r . 
Actually this is not the case; the value U of the internal energy is well defined in the 
entire real axis. Indeed, consider a given parametrization of 7 with variable r G [0, 1] ; 
since the spin connection takes values in the algebra of the Lorentz group, the path- 
ordered integral of the spin connection in the interval [0, r] defines a continuous map 
-0 : T I— > e SO(2, 1) where S0{2, 1) is the infinite cyclic covering of SO{2, 1) 
[10]. By construction, one has ipo = identity and the projection of ipi into S0{2, 1) 
coincides with H{'-f) . The fact that the value of U is well defined in the whole real axis 
should not be surprising; something similar occurs with the time defined by a watch. 
Even if the position of the minute-hand is specified by an angle which is defined modulo 
27r , there is no ambiguity at all in determining, by continuously observing the watch, 
the time interval elapsed with respect to a given moment. 

5. Numerical results. In order to complete the description of the computation 
rules, we need to discuss the specification of the initial data which consists of three 
steps: 

1. it is assumed that all the particles have the same mass m ; a set of values { g^^-^ } 
for the particle momenta is fixed by choosing the spatial components { } 
randomly within a finite domain; 

2. for a given set of values { g^^^ } determined in point 1, the corresponding total 
3-momentum is computed by means of the relation (4.1); 

3. if the total momentum is time-like, by means of a Lorentz transformation A the 
particles system is taken to its rest frame. A is fixed according to equation (4.3) 
which also determines the internal energy U of the system. The set of particle 
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momenta with respect to the system rest frame {p^j^ = ^"'b^'^^ij } represents the 
ordered set of initial particle momenta on which the time evolution acts. 

In the various numerical simulations that we have performed, the number of par- 
ticles has been chosen within 10^ < N < 10^ and the values of the masses belong to 
the range < m < 10^^/46*. Since the total 3-momentum has been determined by 
means of equation (4.1), the values of the mass and of the initial momentum of each 
particle have been chosen to be "small enough" so that no 1/4G ambiguity resulted in 
the computation of the total energy. Instead of performing the Lorentz transformation 
mentioned in point 3, in several occasions we have used a different method to bring the 
system in its rest frame; namely, to the set of the initial particles we have simply 
added a new particle with a suitably chosen momentum so that the resulting spatial 
components P(tot) vanish. This last method is convenient because it avoids strong 
correlations among initial particle momenta. 

5.1. Open universe. We have done simulations consisting of approximatively 4x10^ 
random interactions and we have verified that, as long as U < l/AG , the simulations 
always reached a stationary state within the first 10% of the interactions, which we 
excluded from the final data. The numerical results are expressed by means of the 
equilibrium distribution function /(p) such that 

dN = f{p) d^p 

represents the particle density in momentum space. The normalization condition is 

I d'p f{p) = N . 

Since f{p) is isotropic, it only depends on \p \ ; so it is convenient to introduce a 
distribution function p{p'^) for the energy p^ = ^Jvrfi + \p^ of a single particle: 

P(p°) = ^/(IpI) = ^/(vW^) • 

The p normalization is 

/• + 00 

/ dp^ p{p^) = 1 . 

J m 

When U < 1/46*, we observed that the system always reaches a stationary state. 
Two examples of computed distribution function p are shown in Figure |^ and ^. In 
both cases, = 10^ and m = . Figure |l| refers to the case in which AGU = 
1.296 X 10"^ , whereas in Figure || one has AGU = 9.9968 x 10~^ . In all the figures, 
energies and momenta are expressed in units of l/SvrG. The qualitative behaviour of 
the gravitating gas does not depend on the value of m provided Nm < 1/4G . 
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Figure 1: plot of at N = 1000, m = 0, 4GU = 1.296 x 10"^ 
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Figure 2: plot of at = 1000, m = 0, 4G'C/ = 9.9968 x 10"^ 
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In the low energy limit U <^ l/AG the equilibrium distribution f{p ) is well fitted 
by the Boltzmann distribution that we shall derive later. This is in agreement with 
the fact that, for G — > , the energy conservation law reduces to the classical one. 

As U grows, f{p) differs from the Boltzmann distribution substantially. The high- 
energy tail of p(p°) is rather well fitted by the function a/(p° + b) with a truncation 
at a certain p° = . 

A new phenomenon happens when the value of U approaches 1/4G : the whole gas 
strongly correlates, all the particle momenta tend to become parallel and the overall 
direction fluctuates randomly as the system evolves. Figure |] shows the mean value 
( p ) of p in a finite time interval for AGU almost equal to unity. The time evolution 
presents large fluctuations in |p | which appear as spikes in the picture. 

120 I 1 1 1 1 1 1 1 
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-60 I 1 1 1 1 1 1 1 

-20 20 40 60 80 100 120 

Figure 3: plot of (p ) at = 100, m = 7.161 x IQ-^ AGU = 9.9980 x 10" 
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5.2. Spacelike total momentum. In this case, numerical simulations show that 
the system is unstable and no equilibrium distribution exists. During the evolution, 
the energies of the particles tend to increase overcoming whatever fixed value. The 
way energies diverge is quite peculiar: at first, all energies tend to increase slowly and, 
after a while, they grow faster and faster. The higher the energy is, the quicker it 
increases (positive feedback) so nothing can prevent energies to diverge. At this stage, 
Gott pairs are created and a "jet" appears. Two consecutive particles form a Gott pair 
[11] when the non-Abelian composition p of their momenta, 

exp(p(i)) ■exp(p(2)) = exp(p) , 

is space-like (p)^ = p°'Pa < . A jet is a set of a few consecutive particles -usually 
less than 10- whose momenta are almost parallel and whose energies are much bigger 
than the gas average. Jets have a remarkable property: when particles in the same jet 
interact, their momentum is in general strongly boosted in the jet direction and their 
energies can increase even by 10 -r 30 orders of magnitude in a single transition. Note 
that, as demonstrated in [12,13], Gott pairs cannot be created when AGU < 1 . 

Simple qualitative considerations show that a universe with spacelike total momen- 
tum is either exploding or imploding. Results obtained from our simulations of such 
universes always correspond to the imploding case, where particles looping around 
the whole universe always get boosted in the same direction as their velocity. In an 
exploding universe, particles with speed lower than the global expansion rate cannot 
complete a loop around the universe but, by construction, our algorithm does not take 
into account this global effect. The instability that we find for spatial total momentum 
is in agreement with the big-crunch phenomena studied in [6] . 

5.3. High energies. For a spacetime of topological type 5*^ x IR , the Gauss-Bonnet 
equality implies AGU — 2 . Suppose that one or more particles are very far from all 
the others and do not interact with them; these particles can be neglected so that it 
makes sense to consider a particle system with energy 1 < AGU < 2 . The particular 
case AGU = 1 also corresponds to a (globally static) cylindrical universe. 

In the whole range 1 < AGU < 2 we observed the same behaviour as for spacelike 
total momentum: Gott pairs, jet formation and instability due to divergent energies. 
This suggests that, in our approximation, a closed universe usually develops Gott pairs. 

Because of the finite precision of numerical simulations, it is impossible to constraint 
U to be exactly equal to a predefined value X . Using double precision fioating point 
numbers (8 bytes) we can only obtain simulations with \U — X\ < 10~^^ and U 
fluctuates in that range as the simulation proceeds. We have even seen, though very 
rarely, simulations with AGU < 1 explode due to AGU becoming > 1 by finite 
precision errors. For this reason, it is not possible to distinguish the gas behaviour at 
AGU = 1 (or AGU = 2 ) from the behaviour at very near values of U . What we have 
actually seen is a strong instability in both cases. 
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6. Cleissical statistics. In the next section we shall compute the entropy of the 
gravitating gas when the energy satisfies AGU < 1 and we shall compare the outcome 
of the simulations with the entropy of a non-gravitating relativistic Boltzmann gas in 
(2+1) dimensions. Here we recall a few classical statistical mechanics results. 

For a classical relativistic gas with 2-dimensional volume V , the Maxwell-Boltzmann's 
distribution /o takes the form 



V J 27rkT{m + kT) \ kT 



and has normalization 



/ 



(fx(fph{p) = TV 



The pressure P can be computed in the usual microscopic way: one has to sum the 
variation of momenta over the particles hitting a unitary surface. The result is the 
same as for a perfect gas 

PV = NkT . 
The internal energy turns out to be 

U = U{T,V) = Nip"") = 7VA;rf '"^yj ^ + Nm . 



m + kT 

According to Boltzmann's if -theorem, the entropy is 



S = -kj d'x d'pf ologifo) = kN (^ '^'\_^^^ - log y + log {kT{m + kT))^ 



which satisfies the fundamental relation 

dS{U, V) 1 



(6.1) 



dU T ' 

When m — one obtains 

U = U{T,V) = 2NkT (6.2) 

S S{U, V) ^ kN(^2-\ogy + 2 log . (6.3) 

Note that equation (6.2) is in agreement with the Equipartition Theorem because the 
dispersion relation takes the form ^/\p\'^ ■ 

7. Statistical properties of (2-|-l) dimensional gravitating gas. Numerical 
simulations showed that the open universe with energy iGU < 1 is stable, i.e. the 

system always reaches a stationary state, therefore we shall now examine the statistical 
properties of such a gas. In particular, we will study the entropy and the temperature 
of the system. 
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In order to verify that the stationary state corresponds to thermal equihbrium, we 
have checked that the connected part of the 2-particle probabihty distribution 

f2c{Pi,Pj) = f2{Pi,Pj) - f{Pi) ■ f{Pj) 

is numerically negligible. Wc again apply Boltzmann's //-theorem to give a numerically 
usable definition of entropy; considering only particle momenta, one has 

S = -kj d'pf{p)\ogif{p)) . (7.2) 

For the plots, a "reduced" definition of S has been used: 
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Figure 4: plot of ^ at = 1000, m = 

As one can deduce from Figure 4, the values of S that we found are very well fitted 
by (X log ([/)+ const for energies up to AGU < 0.8 . As U grows, dS/dU decreases 
more slowly than expected. At a certain value U = U , dS/dU stops decreasing and 
S has a flex {d'^S/ dU"^ = 0) . For U > U , the second derivative of S* is > and 
increases as AGU — > 1 . In the same limit, the entropy S seems to diverge. This 
behaviour is supported by the following argument. For a given U , let fm be the 
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maximum value of |) ■ At AGU > 1 we already said that simulations diverge, and 
then it is reasonable to assume that 



lim fm = ■ (7.3) 
Prom the definition (7.2) of S we get 

S - -kj d'pfip) l0g(/(p)) > -k \0g{fm) j d'pfip) = -Nk \0g{fm) , 

then by means of (7.3) one finds 

lim S(U) = +00 . 

Finally, for AGU > 1 there is no equilibrium state, so S is not defined in that range. 

7.1 Temperature The standard interpretation of the entropy behaviour as a function 
of U uses the relations 

1 ^ dS{u,v) 

T dU ^ ' 

and 

1 _ , d-S{U,V) 

C ~ • ^^-^^ 

According to these equations, the temperature T has a maximum at U = U and 
the heat capacity C diverges aX U = U and is negative at U > U . Such results 
look rather strange as they violate the thermodynamical inequalities. In a standard 
statistical system with additive S and C , the violation of these inequalities implies 
instability. This is not true in our system; when AGU < 1 the gravitating gas reaches 
the equilibrium. The standard instability theorems do not apply to our system because 
S and C are not additive. 

In order to illustrate the peculiar properties of the gravitating gas, consider for 
example two identical systems, each one in equilibrium, which are isolated. When 
these two systems are put in thermal contact, in general the resulting system is not in 
equilibrium. 

The temperature T defined in equation (7.4) is not a convenient variable to deduce 
the stability properties of the gravitating system. For the Boltzmann gas, one can give 
different definitions of temperature that are all equivalent. In the case of a gravitating 
gas, these definitions are no longer equivalent. Firstly, we have two energies that 
characterize the equilibrium state: N£ = N{p^) and U. Due to the peculiar energy 
conservation law in (2+1) gravity, NS and U are not trivially related. Secondly, we 
shall introduce the temperature in two ways: by means of appropriate Equipartition 
Theorem (6.2) or by means of the entropy (6.1). For simplicity we shall concentrate 
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on the m — case. By taking into account all the possibilities, the various definitions 
of temperature are: 

T = A T = 

' 2k ' ' 2Nk 



N 



dS_ 



T4 coincides with definition (7.4). 



dS_ 
dU 



kT 




Figure 5: plot of kT at N = 1000, m = 



As we can see from the plot, all the definitions of T are in good agreement with 
each other as long as AGU 0.8 , while at higher energies each definition of T has a 
different behaviour: when plotted as a function of £ , T3 = N{dS/d£)^^ grows more 
than linearly, while 74 = (dS/dU)"^ has the maximum we already anticipated, then 
drops quite quickly. The other two ones Ti = S/2k and T2 — U/2Nk have quite 
obvious behaviours: the first one is linear, while the second asymptotically approaches 
l/8GNk as AGU 1. 

By plotting S as function of logS (see Figure 6), one notices that the fit jS" oc log£^ 
is quite reasonable at all energies. This suggests that T3 is in some sense the most 
suitable definition of temperature. With respect to T — , the AGU 1 limit 
corresponds to S ^ +00 and simultaneously T +00 ; moreover, the heat capacity 
is always finite and positive. 
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Figure 6: plot of ^ at = 1000, m = 



8. Conclusions. A system of gravitating particles in (2+1) dimensions has peculiar 
thermodynamic properties due to the non-Abelian structure of the energy-momentum 
conservation law. Depending on the total 3-momentum, numerical simulations show 
different system behaviours. When the internal energy U satisfies AGU < 1 , an 
equilibrium state is reached; on the other hand, for space-like total momentum or when 
4GC/ > 1 , dynamical instabilities appear. This is consistent with the interpretation of 
particle diverging energies as a big-crunch phenomenon. 

In the presence of thermal equilibrium, we started with the distribution function in 
momentum space and studied the thermodynamic properties of the gravitating parti- 
cles. In the low energy limit, which is equivalent to the low gravity limit, the system 
behaves as a Boltzmann relativistic gas. For 4GC/ — > 1 , the entropy and the mean 
value of the kinetic energy seem to diverge. 

Starting from the various equivalent definitions of temperature in Boltzmann gas, 
we have considered their behaviour also in the energy range of large gravity effects. 
These temperatures turn out to be different. One of them is well behaved also in the 
4GC/ — > 1 limit and does not violate the thermodynamic inequalities; in particular, 
the associated heat capacity is always finite and positive. 
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